Diverging Entanglement Length in Gapped Quantum Spin Systems 
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We prove the existence of gapped quantum Hamiltonians whose ground states exhibit an infinite 
entanglement length, as opposed to their finite correlation length. Using the concept of entanglement 
swapping, the localizable entanglement is calculated exactly for valence bond and finitely correlated 
states, and the existence of the so-called string-order parameter is discussed. We also report on 
evidence that the ground state of an antiferromagnetic chain can be used as a perfect quantum 
channel if local measurements on the individual spins can be implemented. 
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The fields of Condensed Matter and Quantum Infor- 
mation Tlieory share a common interest in the study of 
quantum states of many-body systems. Much of the cur- 
rent effort in Quantum Information Theory is devoted to 
the description and quantification of the entanglement 
contained in quantum states in general: this intriguing 
property of Quantum Mechanics is the basic resource of 
most of the applications in this field, including quan- 
tum communication and computation. Condensed mat- 
ter theory, on the other hand, is deeply interested in the 
strongly correlated states appearing in certain materials 
at very low temperatures, since they describe a variety of 
fascinating phenomena, like the ones occurring in quan- 
tum phase transitions and superconductivity. 

Despite the fact that the number of parameters to de- 
scribe a quantum state scales exponentially in the num- 
ber of particles, it is sometimes possible to capture the 
most relevant physical properties by describing these sys- 
tems in terms of very few parameters. In the case of 
spin chains, for example, two-particle correlations play a 
fundamental role. They allow us to understand several 
complex physical phenomena, like phase transitions and 
the appearance of a length scale in the system. Much in- 
sight has also been obtained by studying exactly solvable 
models such as the AKLT-model which illustrates 
the appearance of the Haldane gap |3] in spin-1 antifer- 
romagnets and the associated finite correlation length. 

From the perspective of entanglement theory, the pres- 
ence of two-particle correlations between two distant par- 
ticles in a many-particle pure state guarantees the pos- 
sibility of establishing EPR-type entanglement between 
them by doing local measurements on the other particles 
0. On the other hand, highly entangled multiparticle 
pure states typically have reduced two-particle density 
operators close to the maximally mixed state and there- 
fore only exhibit very small correlations. Correlation 
functions can therefore only reveal partial information 
about the long-range quantum correlations that ought to 
be present in a state. The so-called Localizable Entan- 
glement (LE) between two particles is defined as the 
maximal possible bipartite entanglement that can be lo- 



calized between them, on average, by optimizing over all 
possible local measurements on the other particles. The 
LE has a very clear operational meaning as it quantifies 
the amount of entanglement that can be localized at e.g. 
the end points of a spin chain and that could be used, 
e.g., as a perfect quantum channel. Just as correlation 
functions induce a correlation length in a lattice, the 
LE induces a new length scale which we call the entangle- 
ment length S_E- It has been proved in jj] that > ^c! 
therefore a diverging correlation length at e.g. a quan- 
tum phase transition implies a diverging entanglement 
length. In the case of ground states of typical spin-1/2 
systems such as the Ising chain and the Heisenberg an- 
tiferromagnet in a magnetic field, we observed that the 
bound is tight = A] ^-^d hence the opposite also 
holds true in this case. This triggered the quest for a 
phase transition that is detected by a diverging but 
for which the S,c remains finite. The natural candidates 
are the spin-1 Hamiltonians that have a Haldane gap 
and hence finite correlation length. 

In this paper, the following results are established: 1/ 
We show that for a family of interesting spin-1 Hamil- 
tonians, including the celebrated AKLT-model, the en- 
tanglement length diverges as physical parameters are 
changed, whereas the correlation length remains finite. 
2/ We calculate the LE exactly for a whole class of 
finitely correlated states Q , which are generalizations of 
the AKLT-ground state. This is a highly nontrivial result 
as the definition of the LE involves a complex optimiza- 
tion over all possible measurement strategies. 3/ The 
so-called string-order parameter 0], reflecting a myste- 
rious topological hidden long-range order in spin-1 anti- 
ferromagnets, is given a natural interpretation in terms 
of the LE. It is shown that, in the family of deformed 
AKLT-modcls, there exist ground states with infinite 
but vanishing string order parameter. 4/ We report on 
numerical results indicating that the LE between the 
two end points of a Heisenberg spin-1 antiferromagnetic 
chain is the maximal possible one. These results indicate 
that ideas and techniques developed in the last few years 
in the field of Quantum Information Theory prove useful 
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FIG. 1: (a) The Valence Bond ground state of the AKLT- 
model each edge represents a singlet, and the dashed cir- 
cle corresponds to a projection onto the symmetric subspace. 
(b) A deformed AKLT model where the singlets are replaced 
by non- maximally entangled states 

to analyze and understand certain aspects of the field of 
Condensed Matter (see also 0)- 

Let us start by recalling the AKLT-model 0, which 
plays a central role in the understanding of gapped spin 
systems. To make calculations simpler, we will con- 
sider an open chain of N spin~l particles at positions 
1, 2, . . . , TV, and with spin-1/2 particles at the ends (lo- 
cations and -I- 1). The AKLT-Hamiltonian is 

AKLT ^ AKLT ^ f " 1 ^ ^ 2 

^AKLT _ ^^AKLT _ ^ ISkSk+1 + -^{SkSk+lf + - 

with S the three spin operators. Each term X^f^J 
is a projector on a 5-D subspace of the 9-D space of 
two spin-1 particles. The ground state can be obtained 
by representing each spin-1 (l...A^) by 2 spin-l/2's 
(1,1,2,2,...) and project them on the 3-D symmetric 
subspace (see Fig. P). It can indeed be checked that 
j^AKLT 1^ orthogonal to all operators of the form 

All <^ A22[Xi <E) |/)i2(/| <E) X2]Al^ ® ^2 (2) 

with Xi,X2 arbitrary 2x2 operators, |/) the singlet of 
two qubits |/) = |01) — |10), and A the 3x4 operator 
that projects a system of two qubits onto its symmetric 
subspace. The ground state of 7jA.klt jg miiq^g ^ ^nd 
can be written as 

1^) = K=iAk) \I)-oi\I)l2 ■ ■ ■ \I)nn+i- (3) 

Indeed, all reduced density operators of two nearest 
neighbor spins of | V) are of the form (|2Jl . 

Due to the special structure of this state, the following 
basic properties of singlets can be used to get some in- 
sight about the quantum correlations present: i/ Given 
a complete spin-1 basis |/3) and a 3 x 4 operator A, then 
there always exist 2x2 operators A^^ such that 

(ft|A=(/|,,Aft®l2. 

ii/ Qubit operators can travel through singlets as 

Al^ (g) =l(g> ayiAl^fayll). 

iii/ The concepts of quantum teleportation and entan- 
glement swapping 8, 9] allow two unentangled particles, 



each of them entangled to a different auxiliary particle, 
to become entangled by doing a joint measurement on 
these auxiliary particles in the Bell basis: 

kk{i\ {\ihk(^\i)kj) =-\ihj 

The crucial observation is now that a von-Neumann 
measurement on the fc'th spin-1 in the basis 

{\ak)} = {|0)fc, \±)k} ^ mk, (I - l)k ± \l)k)/V2} (4) 

exactly corresponds to a Bell measurement on the two 
qubits k and k. If all spin I's are therefore measured in 
this local basis, the mechanism of entanglement swapping 
will produce a Bell state between the two qubits at the 
end points of the chain, independent of its length. This 
proves the existence of an infinite entanglement length in 
ground states of gapped spin Hamiltonians. 

To be more precise, the three properties above allow to 
rewrite \V) in the very convenient matrix product form 
|5||, from which all correlation functions and the LE can 
be calculated exactly. Indeed, inserting a resolution of 
the identity I3N — <8'^i(X]q |Q!i)(ai|) in expression ||2J), 
it follows that 

\V)^ J2 |ai)...|aAr)(l®A""...A"0|/)o,iv+i 

Ql,a2----,C«JV 

(5) 

If one chooses the basis (@J), then A'^ — 0^,^^ = 
(Ty, A^ = with a a the Pauli matrices. The information 
about the remaining Bell state after measuring the spins 
1 . . . A'^ can of course be deduced from the measurement 
outcomes. 

In the case of the AKLT-model, the so-called string 
order parameter has been studied extensively because 
it reveals a hidden topological long-range order. It is de- 
fined as the expectation value of the multi-site observable 

(ao^(®Liexp(z^5fe^))f7^+i). (6) 

The existence of such an order can easily be under- 
stood from the formalism presented here. The operators 
exp(i7rS'^) are all diagonal in the basis l@J. If one mea- 
sures all spin I's in that basis, the fact that the final Bell 
state will be |00) ± |11) or |01) ± |10) is solely determined 
by the parity of the number of times A'q the measure- 
ment outcome is 0; indeed, A"™ • • • will be diagonal 
if and only if A^+ -I- A^- = A^ — Aq is even. The operator 
® exp(z7rS'^) is exactly keeping track of this parity. Tak- 
ing the trace can be interpreted as averaging states after 
measurement, where (8) exp(i7riS'^) introduces a negative 
weight to the states with odd N — Nq. Therefore the 2- 
qubit operator Tri...Ar[((g)^i exp{iTTS^))\V) {V\] has diag- 
onal elements [1/4, —1/4, —1/4, 1/4] and hence maximal 
correlations in the zz-direction, independent of A^. The 
string order parameter is therefore a manifestation of the 
symmetries in the mechanism of entanglement swapping. 
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Let us next investigate what will happen to the en- 
tanglement length when the AKLT-Hamiltonian is de- 
formed. We introduce the following 1-parameter family 
of Hamiltonians: 



N 



(7) 



fc=0 



= ((st)-i®st+i)^w((st)-'®st+i) 

where is defined as 

^ Ife + (cosh(0) - l)Sl + sinh((/))(5^)2. 

The AKLT-model corresponds to = 0, and the per- 
turbation breaks the 0(3) rotational symmetry to 0{2). 
In the limit of — > ±oo, the unique ground state is a 
product state with all individual spins eigenstates of Sf, 
with eigenvalue 0. The unique ground state is completely 
specified by replacing A in (jS)) by 



A = 






exp(— i; 



As will be explained later, one can calculate the corre- 
lation functions and the LE exactly for the end points 
of the chain. The associated correlation length and 
entanglement length S^e are given by 



= 1/ln Vcosh(20)2 -I- 3 cosh(20) 
iE = 1/ln 



^cosh(2(/))2 + 3 + cosh(20) 



and shown in Figure [3 The quantity remains always 
finite and attains its maximum for = 0. is always 
strictly larger than ^c, indicating the presence of two 
different length scales in the ground state. Moreover, 
diverges at the AKLT-point <j) ^ Q. This quantum 
transition with a diverging length scale remains clearly 
undetected by the properties of the correlation functions. 

The optimal measurement basis turns out to be inde- 
pendent of (j) and given by Q. The reason that is fi- 
nite for (/) ^ is the fact that the perturbation effectively 
replaces the singlet |/) = |01) — 110) in equation © with a 
non-maximally entangled state exp((/))|01) — exp(— 0)110) 
(see Fig. nj; using such a state for entanglement swap- 
ping degrades the entanglement (a 6ms is added at each 
step), hence giving rise to an exponential decay and a 
finite entanglement length. 

Let us now develop the mathematical formalism to cal- 
culate correlation functions and the LE. We will consider 
a generalization of the AKLT-states, the family of so- 
called finitely correlated states (FCS) 0, which are, in 
the appropriate limit, dense in the subspace of all transla- 
tional invariant states. These FCS are completely param- 
eterized by a matrix A as in equations H3I5|I . but instead 
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FIG. 2: Correlation length (solid) and of the Entangle- 
ment Length (dotted) for the ground state of the deformed 
AKLT-Hamiltonians H{(j}) O. 



of spin-1 systems we consider general spin-S* systems. 
The spin-l/2's are replaced by a spin-(Z?— l)/2, and |/) 
becomes a maximally entangled state in a D x _D Hilbert 
space [A is a (25* -I- 1) x matrix; note that S and D 
are two independent parameters]. These are all unique 
ground states of gapped local Hamiltonians that can be 
constructed by calculating the orthogonal complement of 
an expression equivalent to (O. For simplicity, we will 
again consider a spin chain with spin S at sites 1 . . . 
and spin {D — l)/2 at the end points. Expectation val- 
ues of the form (t^|Oi ® O2 ® ■ ■ ■On\V) can readily be 
calculated by defining the real x matrices R{0) 
as 



R,.j[6) = M.fcTr {A^6A)aj ® 



(8) 



with {ctq.} a complete orthonormal basis for hermitian 
operators including ctq = 1. The matrix M depends 
on the choice of |/) and of the basis {(Tq}; in the 
case of the singlet state and {tTa} the Pauli matrices, 
M = diag[l, — 1, — 1, — 1]. For the example presented in 
the previous section, e.g., i?(i) is given by 



i?(l) 



/ 3cosh(20) 
-10 
0-1 
\-3sinh(20) 



sinh(20) \ 



-cosh(2(/)) j 



Evaluating expectation values becomes equivalent to cal- 
culating the (0, 0) element of the matrix product 



i?(Ojv)-R(Ojv-i)---i?(Oi 



(0,0) 



The normalization of \y) is clearly given by 
[(i?(l))^]^pg^, and therefore correlation functions 
can be calculated as 



(Oi ® 62 



■■■On) - 



i?(Ow)---i?(Oi) 



(0,0) 



mt)r\ 



(0,0) 



To calculate the LE, we need to do an optimization over 
all possible local bases (including overcomplete ones 
corresponding to generalized measurements). If one mea- 
sures all spins k = 1 ... iV in the basis {|/3fe)}, the state 
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of the two extremal spins is conditioned on the outcomes 
. . . , Pn} and given by 



(9) 



where we used the notation {Pi\A = {I\A^' ® 1. The 
probabiUty for this outcome to happen is 



(V'ft,...,/3wlV'/3i,...,;3») 
{V\V) 



FoUowing ^3 J a generahzation of the concurrence 
for pure bipartite D x D states is given by C(|x) = -B (8) 
1d|/)) = \det{B)\'^^° /{x\x)- Using this measure, the 
average entanglement factorizes and is given by 

Ep^^^JdetiA^- ■■■AP^)\i _ n. (Eft|det(A'3>)|* 



{V\V) 



{V\V) 



Obviously, the optimal strategy is to measure all spins 
in the same basis that maximizes J^p- I det(A'^')|"D . This 
problem is equivalent to calculating the entanglement of 
assistance (EoA) of the (unnormalized) x state 
A''' A, and can in general easily be done numerically. 

In the case of singlet valence bonds {D — 2) and A an 
arbitrary 3x4 matrix, this EoA can be calculated exactly 
[l3| . Making use of the results in one obtains the 
exact expression for the LE: 



^\^,^{MR{l)M{Ril))T) 



N 



'0,N+1 



mm 



N] 



(10) 



J 00 



Here was defined in I©, M = diag[l, -1, -1, -1], 

and Amax means the largest eigenvalue. 

Let us investigate the general conditions under which 
the entanglement length as defined by the LE H1U|I di- 
verges. It can be shown P| that this will happen if and 
only if there exists an operator Q such that the EoA of 
(Q^ ® Q'')A'' A{Q (g) Q) is maximal, i.e. when it can be 
written as a convex combination of maximally entangled 
states (such as for the AKLT); this is indeed the necessary 
and sufficient condition for perfect entanglement swap- 
ping to be possible. Let us compare this with the condi- 
tion for the presence of a generalized string order param- 
eter, which we define as in © but with an extra opti- 
mization over the hermitian operator X in 0^ exp(i7rArfc) 
instead of Sz- The condition is now the existence of an 
operator Q such that (Q'l' (g) Q'')A'' A{Q (g) Q) is diagonal 
in a basis of maximally entangled states that have either 
support in the |01) ± |10) or in the |00) ± subspace; 
this gives two extra nontrivial constraints as compared 
to the condition for a diverging entanglement length, and 
hence there exist ground states with vanishing string or- 
der parameter, but infinite entanglement length. 

Let us conclude by discussing some generalizations to 
the present results. First of all, the formalism developed 



in this paper allows us to calculate the localizable en- 
tanglement for ground states of arbitrary Hamiltonians 
numerically by making use of the Density Renormaliza- 
tion Group formalism (DMRG) the fixed point of 

the DMRG algorithm yields a FCS 0, on which we 
can apply the machinery presented. Secondly, the results 
presented equally apply to higher dimensional AKLT- 
models, and reveal an intriguing connection with quan- 
tum computation |l7| |. Finally, we have done numeri- 
cal diagonalizations of the Heisenberg antiferromagnetic 
spin-1 Hamiltonian H = SkSk+i with spin 1/2's at 
the end points Q . The localizable entanglement between 
the end points is again 1 independent of the number of 
spins N and is obtained by measuring in the optimal 
basis (@J for the AKLT: the entanglement length in the 
Heisenberg spin-1 antiferromagnet is also infinite, and 
hence this ground state could be used as a perfect quan- 
tum channel. It is interesting to note that the entan- 
glement length for the spin 1/2 antiferromagnet is also 
infinite the presence of a Haldane-gap severely af- 
fects the correlation length, but does not seem to affect 
the entanglement length. 
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